Small Energy Scale for Mixed- Valent Uranium Materials 
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We investigate a two-channel Anderson impurity model with a 5/^ magnetic and a 5/^ quadrupolar 
ground doublet, and a S/'^ excited triplet. Using the numerical renormalization group method, we 
find a crossover to a non- Fermi liquid state below a temperature T* varying as the 5/^ triplet-doublet 
splitting to the 7/2 power. To within numerical accuracy, the non-linear magnetic susceptibility 
and the 5/^ contribution to the linear susceptibility are given by universal one-parameter scaling 
functions. These results may explain UBeia as mixed valent with a small crossover scale T* . 



PACS numbers; 75.20.Hr, 71.10.Hf, 71.27.+a, 72.15.Qm 



The possibility of non-Fermi liquid (NFL) behavior in 
Ce- and U-based alloys has been discussed intensively 
since the discovery of the anomalous temperature de- 
pendence of their resistivity, magnetic susceptibility and 
specific heat coefficient. A candidate for the explana- 
tion of the NFL state is the multichannel Kondo model, 
widely used for the /-shell materials |l). For the /-shell 
impurities, the intra-atomic interactions such as Hund's 
and spin-orbit couplings have to be taken into account in 
the presence of the crystalline-electric field (CEF) , giving 
rise to various types of scattering of conduction electrons 
|g,||. Experimental investigations of dilute Ce or U al- 
loys suggested that single-site effects of such a magnetic 
ion are important for NFL physics. In the U case, the 
observation of a logarithmically divergent specific heat 
coefficient [Q strongly supported the quadrupolar (two- 
channel) Kondo scenario H . As evidence of single-site U 
effects, such a logarithmic anomaly was observed in the 
dilute U limit of Ua;Thi_a;Ru2Si2 (6J, although the resis- 
tivity in this metal cannot be explained by the simple 
two-channel Kondo model 0. 

According to photoemission and inverse photoemis- 
sion studies on U compounds [0 , well-separated atomic- 
like peaks cannot be seen in 5/-spectra, while for Ce 
compounds the peaks can be described by an Ander- 
son impurity model. This implies that 5/-properties 
look much more like mixed valence. Recently, the two- 
channel Anderson impurity model in the mixed-valent 
regime was proposed to account for the NFL physics 
of UBei3 0. In this model, a low- lying Fa quadrupo- 
lar (non-Kramers) doublet in the 5/^ configuration and 
Fg magnetic (Kramers) doublet in the 5/^ configuration 
are taken into account and excited CEF states are ne- 
glected. Both states mix with each other via the hy- 
bridization between the localized /-orbital and conduc- 
tion band. The model successfully describes the temper- 
ature dependence of the non-linear magnetic susceptibil- 
ity observed in UBeia Q and in Ui-^^Th^^Beis |jll|], sug- 
gesting that strong quantum fluctuations should drive the 
CEF state of U ions to a mixed-valent state between U+^ 
and U+"' [Ol. However, except for virtual effects like the 
Van Vleck magnetic susceptibility, excited 5/^ CEF lev- 



els were neglected in the model. This leaves a large ques- 
tion, however: how can the small energy scale (~ lOK) 
apparent in thermodynamic data for, e.g., UBcia be rec- 
onciled with such a mixed valence model for which the 
smallest plausible value is large (~ 150K)? 

In this paper, we show that the inclusion of the excited 
CEF levels resolves all the above controversies. The dy- 
namics of the excited CEF states reduces the crossover 
temperature to the two-channel Kondo state significantly 
even in the mixed-valent regime. Using Wilson's numer- 
ical renormalization Group (NRG) method [Ijjl^l for a 
cubic CEF case, we show clearly that a small energy scale 
T* is related to the competition between different types of 
NFL fixed points, which is very important to explain the 
observed non-linear susceptibility within the two-channel 
Kondo scenario. The crossover in the Kondo effect is vi- 
sualized by flows of NRG energy levels and T* can be es- 
timated by observing the recovery of the symmetry asso- 
ciated with the two-channel Kondo model. We find that 
T* obeys a power-law with respect to the CEF splitting 
and can be much smaller as the splitting decreases. In 
addition, the non-linear magnetic susceptibility is domi- 
nated by the magnetic moment of the excited configura- 
tion, and obeys a one-parameter scaling law set by T* . 
We expect this result for the small energy scale to be 
generic for ions such as Pr, Tm and U which fluctuate 
between configurations with internal degrees of freedom; 
the exponent of the power law will vary with the details 
of the CEF spectrum of the ground configuration. 

In the present work, we study the two-channel Ander- 
son impurity model in a cubic system, taking into account 
the first excited r4 triplet state explicitly as well as the 
lowest-lying V^ quadrupolar doublet in 5/^: both states 
are coupled to the Fy magnetic doublet (the lowest pos- 
sible state in 5/^), via the hybridization V with the Fg 
conduction electrons. Note that the choice of hf^ in our 
work is motivated by calculational simplicity; there is no 
qualitative difference with respect to the 5/^-5/^ picture 
used in Ref . . The corresponding Hamiltonian is given 
by 



H = ^ekcl^ckm +^Em\M){M\ 
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where a single- /-electron state |a) represents the T-; 
states; a two- /-electron state \M) represents the T^ and 
r4 states. The operator c\.^ (ckm) creates (annihilates) 
a Fg conduction electron with a wave number k, the ki- 
netic energy of which is given by Sk ■ The energy of the 
/-electron state Em {— E-p^ , E-p^ ) is measured from that 
of the Tj magnetic doublet. The Clebsch-Gordan coeffi- 
cient Ca,m;M is givcn by (q;|/„i|M), where fm annihilates 
an /-electron with Fg symmetry. If the energy of the F4 
triplet state is so large that the contribution from the 
triplet can be neglected, the Hamiltonian (|^) becomes 
simpler as presented in Ref. f^. Assuming further that 
TTpV"^ <C I^Fsl (p is the conduction-electron density of 
states at the Fermi level), the Hamiltonian is reduced to 
the two-channel Kondo Hamiltonian with 5" = 1/2 {S is 
the size of a local spin) g,^, which has been solved by 
various methods ||l[ . We can label the four Fg states by a 
tensor product of two quadrupolar and two magnetic in- 
dices. Each of them has SU(2) symmetry in the restricted 
Hilbert space of this isotropic Hamiltonian. In our model, 
on the other hand, we have to treat the cubic symmetry 
explicitly because mixing of F3 and F4 states prevents 
the SU(2) labels from being good quantum numbers. 

In spite of such a complicated atomic structure of the 
impurity, the Hamiltonian Im) can be solved by the NRG 
method |13,0]. For the purpose of numerical calculation, 
it is transformed to the following hopping Hamiltonian 
with the recursion relation 

Hn+1 = A^/^i/jv + E^^Wl.m^Wm + h-C), 
m 

Ho = A-'^'[Y,Em\M){M\ 
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Here Hq describes the impurity site, and sjj„j creates a 
Fg electron in the conduction band discretized logarith- 
mically as D/A (D is a half width of the band). The en- 
ergy of two /-electron states Em and the hybridization 
F are equal to Em / D and pV^ /D, respectively, except 
for some A-dependent scaling factors |1^. We diagonal- 
ize the NRG Hamiltonian H^ at each NRG step to solve 
the eigenenergies, and we construct H^+i from TJ^r, us- 
ing the Clebsch-Gordan method for the cubic symmetry 
|l4[ , which speeds up the calculation ten-fold. Introduc- 
ing an axial magnetic field hzfl on the impurity site, we 
can calculate the temperature T dependence of the im- 
purity magnetization by using 

TrJ,exp[-^(ffjv + i/;r)] 



M, 



Trexp[-^(i7Ar+i/;,)] 



(3) 



Here the axial component Jz of the total angular mo- 
mentum is coupled to the magnetic field in the Zeeman 



term H'^ = —gjJzhz {gj is a Lande's g factor and we 
take fiB — kB = 1). The physical temperature T and the 
rescaled magnetic field hz are given by 



h^ = 



1 + A-i 



A(^-i)/2/, 



z,0, 
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respectively, where we take /3 ~ 2. The linear susceptibil- 
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ity xlmp ^^*i non-linear susceptibility xlmp are calculated 

by expanding Mz in hzfl (M = x''^^^ + X^^^h^l^ H )■ 

In the NRG calculation, these are explicitly given by 
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This calculation of the impurity susceptibility gives suf- 
ficient accuracy over a wide range of temperatures, as 
long as the hybridization V is restricted to upV^ <^ D. 
Throughout the NRG calculations, we take A = 3 and 
keep 500 or 800 lowest-lying states at each NRG step. 
Although the latter number reduces the numerical error 
relative to the former, a qualitative difference does not 
appear. We show the former results in this paper. 

First we shortly discuss the case Er^ ^ Ep^ in the 
mixed- valent regime H. As the temperature decreases, 
two characteristic temperatures are obtained: at a higher 
temperature, the highest of the F3 quadrupolar and F7 
magnetic doublets is screened by the conduction elec- 
trons; at a lower temperature, the other remaining dou- 
blet is screened. The NFL behavior due to the two- 
channel Kondo effect appears at the lower temperature 
Tk- Tk increases with the increase of V and the decrease 
of I -Era I- According to the analysis for the realistic pa- 
rameter region of the Hamiltonian, the evaluated Tk is 
larger by an order of magnitude than the lOK observed 
in UBei3. If the contribution from the F4 triplet to the 
Kondo effect is taken into account, Tk is enhanced as the 
CEF splitting A = Ep^ — Ep^ (> 0) approaches zero, and 
a smaller temperature T* appears at which we find the 
crossover from the CEF states to the NFL states. Be- 
low we discuss this energy scale, restricting ourselves to 
the realistic case where /^ is the most stable configura- 
tion {Et^ < Er^ < 0) and the two-channel Kondo effect 
occurs at low temperatures. 

Our NRG results show that the stable fixed points are 
classified into two types: for the first, the low-lying exci- 
tations can be described by assuming a F3 local moment 
coupled to the conduction electrons. For the second, a 
r4 local moment is dominant in the low-lying excitations. 
The stability of the fixed points depends on whether the 
bare CEF splitting A is above or below a critical value Ac 
(> 0). The following discussion is restricted to the case 
A' = A — Ac > where we obtain the first fixed point. 
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FIG. 1. Flow of the most relevant NRG energy levels 
for the odd number of NRG steps A*' and A = 3. We take 
f = 0.15, £;r3 = -0.200 and Er^ = -0.168 here. Each 
energy level is labeled by {Q,ri), where Q represents the to- 
tal number of particles measured with respect to the ground 
state. At the fixed point, the ( — IjT^) and (IjFe) doublets 
have a center of gravity energy of 0.2, while the degenerate 
(0, r4) and (0, Fs) triplets have an energy of 0.8. These val- 
ues correspond to the first and second excited NRG energies 
given by the two-channel Kondo model for A = 3 [ISl. This 
implies that only particle-hole symmetry breaking causes the 
splitting of the (—1, F7) and (1, Fe) doublets and that of the 
(-2,F3) and (2,F3) doublets. 



described by the two-channel Kondo model, as is easily 
found from a flow diagram of the low-lying NRG energy 
levels in Fig. 1. For a small number of NRG steps, the 
almost degenerate (0,F3) and (0,r4) states correspond to 
the lowest-lying CEF states with small splitting A. As 
the number of NRG steps increases, the splitting of (OjFa) 
and (0,r4) levels becomes larger, and the former becomes 
the lowest state at the fixed point (where all the NRG 
energy levels reach constant values). The double degen- 
eracy of F3 corresponds to that of a local spin of the 
two-channel Kondo model M. Practically we define T* 
as the temperature corresponding to the number of NRG 
steps at which the difference of the NRG energy levels 
between the (0,F4) and (OjFs) triplets reduces to 0.01. 
The degeneracy of these states is one piece of evidence 
for the two-channel Kondo fixed point. 

Figure 2 shows that T* is proportional to the power 
law with an effective CEF splitting A*. In analogy with 
Boltzmann weight factors, A* is associated with the oc- 
cupancy weights of /-electrons n^ and n1 in the F3 dou- 
blet and F4 triplet at the fixed point, respectively: 



A* = An In 



3nq 
2^ 



(6) 



Here Aq depends on only the hybridization F if the other 
parameters are fixed. When A* goes to zero with A', 
both the F3 doublet and F4 triplet become relevant with 
the same weight. Each line in Fig. 2 has almost the same 
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FIG. 2. Relation between the small energy scale T* 
and the effective CEF splitting A* for f = 0.10 (o), 0.15 
(D) and 0.20 (O), where we fix Ev^ at -0.200. The corre- 
sponding lines are represented by T* = a(A*/Ao)'', given as 
(a,6) = (2.7xl0~^3.5), (4.2x 10~^ 3.7) and (5.6xl0~^ 3.7), 
respectively. Inset: A*/Ao versus A' (= 3[A - Ac]/2). For 
the above ordering of F values, 3Ac/2 is 0.018, 0.031 and 
0.044. 



power, implying T* ex (A*)^/^. As A' approaches zero 
(we use the NRG scale A' = 2/[l + A"i]A' in Fig. 2), it 
tends to be proportional to A* and satisfies T* ex (A')^/^. 
This relation can be derived from the equation for a small 
energy scale Tggp ~ ro(A/To)^/4 appearing in Sec 5.3.2 
of Ref. PJ . Since Tq is the temperature at which the F3 
doublet is screened by the conduction electrons, written 
as To = (D/A)3/2TK,r3, where Tk^F:, is the Kondo tem- 
perature for the case without the higher F4 triplet [[l5[ , 
we obtain 



-'CEF 



1 D^ 



4Ti 



K,r3 



7/2 



(7) 



This resembles the relation of T* and A' if A' is small. 
Actually the hybridization renormalizes A to A' as ob- 
tained by our NRG calculation. For the two-channel 
Kondo effect, TR.ra is given by pV^^a^^iEY.^jlpV'^^, 
where the ratio V^g/V'^ = 8/21 comes from a Clebsch- 
Gordan coefficient when we derive the effective exchange 
interaction due to the F3 doublet. Using Eq. (M), we can 
roughly estimate the ratio T*/{A'y/'^ as l/(lG^TK,r3), 
where D is taken to be unity. 

In Fig. 3, we show the temperature dependent mag- 
netic susceptibility rescaled by T* . The linear suscepti- 
bility Ximp i^ approximately given as a sum of Zeeman 
and Van Vleck parts. The former term is dominated 
by the F7 magnetic susceptibility Xi 



XI 



(1) 

impT? ' 



Each curve 



Ximp.r. is wefl fitted by a function f{T/T,)lT,, loga- 
rithmic below T ~ Ts (Ts ~ lOT*). On the contrary, 
such a simple scaling function cannot be found for the 
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FIG. 3. Temperature dependent impurity magnetic sus- 
ceptibility. The curves are for F = 0.10, Er^ — —0.200, and 
Br4 = -0.160 (O), -0.165 (+), -0.170 (□) and -0.175 (x), 
for which T^ (~ lOT*) is, respectively, 3.6 x 10-^ 2.3 x 10-^ 
1.1 X 10~^ and 1.7 x 10~*, and the 5/^ occupancy in the 
ground state is 0.775, 0.765, 0.751 and 0.727, respectively. 
Upper curve: The 5/^ contribution to the linear susceptibil- 
ity Ximp- Lower curve: The total non-linear susceptibility. 



Van Vleck contribution, which is expected to be given by 
a — by^T/Ts 0]. In this case, a and b appear to be com- 
plicated functions of T* and A*. Since JiFal J^|r4)|VA 
is very large in our model, mainly the Van Vleck part 
contributes to Ximp' ^'^^ ^^'^^ ^ simple scaling relation 
for Ximp cannot be found. 

On the other hand, the curves for the non-linear sus- 
ceptibility T^Ximp "^o scale with T/Tg, implying that the 
main contribution comes from the F7 magnetic doublet 
even when /^ is more stable than /^ (0.5 < 713 -f n| < 



0.8). Usually the r3 quadrupolar moment gives rise to a 
positive log(TK,r3/7") dependence in Ximp due to the two- 
channel Kondo effect, which should be caused by the Van 
Vleck coupling of the F3 doublet and F4 triplet. How- 
ever, the magnetic part derived from the Fy behaves as 

— 1/T^ -^ — 1/T below Tg and overcome the quadrupo- 

(3) 
lar contribution. This trend of negative Ximp ^^ i^ good 

agreement with the non-linear susceptibility derived by 
the simple model where the dynamics of the F4 triplet 
was neglected [g[. However, for comparison of the scal- 
ing behavior, Tk must be replaced by Tg. For Er^ < 0, 
Tk is approximately equal to Tkxs in Eq. (|j) and is very 
large in the strongly mixed-valent regime. On the other 
hand, T* (~ O.lTg) can be much smaller when A is close 
to Ac. We conclude that T*, characteristic to the com- 
petition between the two different types of NFL fixed 



points, can represent the small energy scale of around 
lOK for UBei3. 

As discussed above, the low-energy physics below T* 
is described very well by the two-channel Kondo model 
as long as the CEF splitting satisfies A > Ac. In the 
case A < Ac, low-lying excitations dominated by the 
relevant triplet at low temperatures are very different 
from those given by the two-channel Kondo model. The 
study of this case is also required and motivated by a 
recent experiment of Ua;Yi_2;Pd3 where a F5 triplet is 
almost degenerate with a F3 doublet Q . We anticipate 
that a small crossover scale will still result for A < Ac, 
albeit with a different power law exponent. Indeed, this 
result should be generic for mixed-valent ions such as Pr 
and Tm, having configurations with internal degrees of 
freedom for the lowest CEF levels. 
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